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SUMMARY \ 

The boundary conditions corresponding to the design problem when the blades being simulated by the 
bound vorticity distribution are presented. The 3D flow is analyzed by the two steps 52 - 51 approach. In 
the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an 
axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of 
the blade producing the flow channel striction is taken into account by the modification of metric tensor in 
the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is 
satisfied automatically. The governing equation is deduced from the relation between the azimuthal component 
of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided 
by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream 
sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is 
taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet 
is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential 
to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this 
condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface 
obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the 
penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the 
inverse problem. The detection of this flux leads to the rectification of the geometry of the blades. 
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Superscripts 


angular velocity of the rotor 
pressure 
total pressure 

stagnation enthalpy or ptfp 
rothalpy or H + w(V$) 
blade force 
dissipative force 
efficiency 

transformed coordinates system in Si approach 

leading edge 
trailing edge 
reference 

tangentiel component 
normal component 
nodal point indices 

upper side of the blade 
lower side of the blade 


1. INTRODUCTION. 

Most of the blading design procedures consider the velocity distribution on both sides of the blade as the 
initial data, the inverse problem becomes ill-posed and the designer loses the control of thickness distribution of 
the blade. To overcome this deficiency, this paper suggests an inverse method by representing the blades by a 
distribution of bound and free vortices which produce the desired swirl (V^r) variation. By introduction of the 
notion of associated elements on both sides of the blade in respect of the thickness distribution, and by imposing 
a conservative flux penetration through each pair of the associated elements when the geometry of the blade 
is not yet well defined, we obtain the well-posedness of the inverse problem. The iterative rectification of the 
camber surface in order to cancel the flow penetration leads to the final geometry of the blade. Treating first 
the 2D cascade design, §2 is devoted to show how to get the well posed inverse problem with the appropriate 
boundary conditions applied on the presumed blade contour, and the procedure leading to the rectification of 
the camber lin e related to the penetrating flux of the fluid determined on both sides of the blade. To treat the 
quasi 3D design, the S2 and Si approach as proposed by C.H. Wu [l] is adopted. The loading produced by the 
velocity difference between the two faces of the blade is directly related to the bound vorticity distribution that 
the blade has to generate. Assuming the number of blades infinite, the vortex distribution as well as the flow 
field become axisymmetric (S2 flow), §3 shows how the blade thickness distribution and the loading distribution 
can be taken into account in this scheme, and how to deduce the pressure distribution on the blades when their 
number is finite. An application to the case of the centrifugal impeller is presented. The loss scheme by the 
introduction of a plausible value of efficiency rj for each streamline as suggested by J.H. Horlock [2] is used. This 
approach opens up possibilities for the elaboration of a design which maintains the assigned value of the total 
pressure gain in each stage by modifying the (Vj^r) distribution in free space between blade rows. §4 is devoted 
to the blade to blade flow (Si) inverse problem, the boundary conditions for 2D inverse problem are transposed 
to this quasi-3D flow. The stream function is used to define the flow field and the finite volume method is used 
to solve the problem. Examples show the results concerning the design of centrifugal impeller. 


2. INVERSE PROBLEM FOR THE 2D CASCADE. 

Figure 1 shows the geometry of the blade characterized by its thickness distribution and the shape of 
its camber line. The arc elements taken respectively on the upper side and the lower side tangential to two 
inscribed circles centered on the camber line at x — dx/2 and x + dx/2 are called associated to the camber line 
element. The center of these associated elements are characterized by the abscissa x of the camber line element. 
Let V Q represent the upstream velocity, h the pitch of the cascade, a and /? the inlet and outlet flow angles, the 
circulation T of the bound vortex generated by the blade is given by: 



r = V G h(sin a — cos a tan ft) 

Tb. bo^ v„ rt « dMbrt. on *. bUd, c,» b. 

monotonic increasing function of x for the inverse problem. / **.) “ the traLling 
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on both sides of the blade [3]: 
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be admitted., In order that the boundary condition does not produce extra , V 

each pair of associated elements is to be conservative, this implies: 
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The solution of the inverse problem determines the flux penetrating through the associated boundary elements, 
the camber line inclination correction Sd is given by. 


6# 


= 0.5[tan 1 ( _ ^) + + tan 1 ( y ( ) ] 


(2.3) 
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S. MERIDIONAL FLOW, S 2 APPROACH. 

, ... first 8teD the vort ex distribution is transformed into an axisymmetric one by spreading it in the 

Let7 x e = 6 > a nd ( 3 represent the body fitted curvilinear coordinates (Fig. 4), the meridional velocity is 
represented by:’ V = V 1 ?! + V 3 «s = IF 1 ft + IF*?,, the continuity equation becomes: 
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where g represents the determinant of the modified metric tensor due to the flow channel striction pr^duced by 
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y to represent tbe flow field by imposing: 
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( 3 . 2 ) 


yfpgdi* " V y/P9 d V 

tbe equation (3.1) is satisfied automatically. Tbe governing equation for ^ i* obtained by writing V x V = n J e 2 , 
represent* tbe axunutbal component of V x V^, it ia deduced from tbe bub to sbroud equ il ibrium 
condition. Let 
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Tbe momentum equation ia: 


Pt 

and 

II 

k 

1 

^ = £f + w(V,r) 

p 


P 

2 

2 


r-v# 


( stator 


\-v; 

P 

P 

\ rotor 


( 3 . 3 ) 


In fact, there ia a preaaure gradient in tbe azimuthal direction in tbe flow apace between blade*, in tbe axial 
eymme’tric S2 flow where tbe number of blade* ia auppoaed infinite, thi* preaaure gradient diaappear* and tbe 
volume force Ft,/p due to tbe blade* ha* to be added in the momentum equation. Tbe lot* acbeme [2] related 
to tbe plauaible value of efficiency rj for each atreamline of the atage ia added, thia acbeme auggeat* that tbe 
diaaipative force F^/p ia related to tbe variation de Vjr via tj : 


f(»J- ■*&*)]? >Utor 

P ~ \ (1- h)i#yrl^ • V(V # r)l# rotor 


( 3 . 4 ) 


— o aa well a* fj = 0 are impoaed in the free apace. Figure 5 ehowa the relation between the kinetic moment 
diatribution in tbe blade row apace and the circulation of the bound vorticea produced by the blade*. Let T^, 
denote the circulation generated by the blade in the section cut by an axiaymmetric stream surface \j> — cte, the 
kinetic moment (Vjr) m ,k generated by the bound vortice* located between the leading edge and tbe abscissa m 
can be represented by: 

(V.r)».* = (V,r)„.* + ^r*/(m,V) (3-5) 

Using (3.4) and adopting that dl/d£ 2 or dH/d( 7 being equal to — ( ) 2 / P the dissipative scheme, the 

azimuthal component of the momentum equation leads to: 

{ f b ) 2 _ r v i B(V$r) 3 9(M 1 ( 3 . 6 ) 

~- [V ee + ae J 


where W'a = V#r + ur 3 and V 3 = Vjr. The coordinates «y»tem C “ choaen ao that the constant £ 5 line* are 
iteratively replaced by tbe streamline*. The component following t s of the momentum equation represent* tbe 
hub to shroud equilibrium condition, which gives: 


1 + (*»)» (£}», 
+ P P J 


Let n design the normal of the camber surface of the blade, we have: 


niT 3 -f n 3 r 2 4 - n s r s 


As W 1 n in the rotor and V 1 n in the stator, we have: 
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U»ing (3.4), (3.6) and the la*t 3 relation*, (3.7) become*: 
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In the free apace, the component following F s of the momentum equation lead* directly to: 


free apace y/gCl 7 = yT { 
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The dot product of the momentum equation with V in the atator and in the free apace or with in the rotor 
lead* to the following relation* which aerve to update the nodal value* of H or I: 
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where 3( )/dm denote* the meridional atrearnwiae tangential derivative. Writing V xU - 0 3 e 2 , we obtain 
the governing equation of yp: 


a (SiL -l g ( j» ) 

dt^py/gdt*’ 


a /inilt _ = /jn 2 


(3.10) 


For the inveree problem, the di*tribution of V t r i* assigned, uaing (3.8), f) i> updated iteratively. Let l *> e 
camber aurface of the blade be defined by t = * a U l ,C 9 ) + cte, if the coordinate line. ^ updated to 

the atreamline* iteratively, ( 3 can be computed uaing the alip condition: 


w« 


(3.11) 


Figure 5 ahow. the geometry of the blading of a multiatage turbopump obtained by wiving the inverae problem. 
The CPU time on IBM 3090 in acalar mode i> about 1 minute for the entire turbopump. The gnd u*ed for the 
S2 computation i* 300x16. Figure 6 ahow* the compariaon of the centrifugal impeller. deaigned with tj « 1 and 
t) < l having the ume level of total preature gain. 

Blade aurface preaaure evaluation. - Uaually the S2 approach lead, to the determination of the mean 
velocity on both facet of the blade: 


( rotor 
\ itator 
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Let AI7 denote the difference of the ab.olute velocitie. (V* - V~) or the relative velocity (*' - * ) °® tbe 

two face, of the blade, when the number of blade, i. finite, thi* difference i. related to the local density of 
bound vortex generated by the blade. In the S2 .cheme, con.ider the blade wet, on cut by a J - cte .urfac* 
the flux of bound vortice. generated by the element 6? of the blade i. determined by the fluoffi through 
Z elementary aurface («)* if = y&Wf. where ^ 3 .hould be equal to LJN>. U.ing the Stoke* region 
that implie* the circulation produced by A V i* equal to the flux of the bound vortice. we get the following 

relation: 
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(3.13) 
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where f) denotes the local angle of the blade section with respect to the meridional plane. (3.12) and (3.13) are 
used to compute surface velocity on both faces of the blades, then the pressure distribution by the S2 approach 
can be deduced. (Fig. 9) 

4. BLADE TO BLADE FLOW, SI APPROACH. 

The blade to blade flow confined in each axisymmetric stream sheet is analyzed in order to define the final 
geometry for each section of the blade and to obtain the pressure distribution. At the beginning, the contour 
of the blade is created from the camber line obtained from the S2 step with the assigned thickness distribution. 
The conformal mapping (m,8) => (x 1 ,! 2 ): 



(4.1) 


transforms the blade to blade flow confined in an axisymmetric stream sheet into a 2D cascade flow in the 
(i 1 ,! 2 ) plane. The body fitted coordinate system constituted by the equipotential lines f 1 = cte and the 
streamlines £ 2 = cte of a Active 2D flow around the cascade is created using the panel method [4j. In this 
system, the continuity equation becomes: 

jjlrt + <«> 


where U x represent the contravariant components of the absolute velocity V for the stator and relative velocity 
W for the rotor and 

_ D(x 1 ,x 2 ) , r ^2 


where D(x l , x 2 )/D(£ 1 , £ 2 ) denotes the Jacobian, r represents the local thickness of the stream sheet. Introduc- 
ing the stream function t b with 

ri_ 1 


U 1 = 


U 2 = - 


Py/gdp 

1 dtp 


(4.3) 




(4.2) is satisfied. Prom the momentum equation, we can show that the free vortex of the absolute velocity 
shedding from the preceeding blade row must be tangential to the axisymetric stream sheet, the governing 
equation of the blade to blade flow stream function is deduced from this condition: for the relative flow around 
the blades of the rotor, we have: 


r d / gw d\p . d , gn dj>_ _ 
l dt l[ Py/gdt' )+ d( 2K p y /gd? ) l 

dg 21 W l , dg 17 W 2 , „ log r 

- ~^T~ + -gp- +2 v^v^r 


Boundary conditions for the inverse problem: 


Flux conservation: 

Bound vorticity assigned: 


= 0 

[VMf 1 -ur^dd]* = rdf 


(4.5) 


The solution of the inverse problem leads to the determination of flux penetration on the blade contour, the 
camber line inclination correction is given by: 




r WV 


(4.6) 
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Figure 7 shows the network (f l , f 2 ) around a blade row for an impeller. Figure 8 shows the comparison 
of the camber lines of the impeller obtained from the S2 approach and rectified by the SI approach. Figure 9 
shows the pressure distributions obtained from the S2 and Si approaches. For the case of the turbopump the 
loading is optimised to avoid the cavitation. The results from the S2 and SI computations are similar, but not 
identical the need of the SI computation to obtain the final geometry definition of the blades is confirmed, or 
one stream sheet, the CPU time on a IBM workstation RISC 6000/320 is about 40 minutes, or about 5 minutes 
on IBM 3090 in scalar mode. The grid used is 150x16. 


5. CONCLUSION. 

The representation of the blades by the vortex distribution enables the formulation of the well-posed 
inverse problem, and which leads to design the blading of a turbomachine. The two steps S2 - Si quasi-3D 
approach has been applied on different axial and radial geometries. Several kinds of loading function have 
been tried The results show that the success of the blading design depends greatly on the meridional ( »r) 
distribution assignment associated with the loss distribution. To optimise the design m order to avoid the 
formation of the cavitation or the separation of the boundary layer in the design condition, when the loading 
is not too high, experiences show that an adequate modification of the bound vortices distribution function / 
may effectively lead to prevent the surface pressure to be lower than the cavitation level or to maintain t e 
adverse pressure gradient below the boundary layer separation criterion. The inverse problem procedure has 
been elaborated to calculate the turbomachines in incompressible range, the research works are planning o 
extend this method to make the transonic designs. 
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Fig. 4. Boby fitted coordinate system £ , £ 




Meridional section of a multistage turbopump. 









The blading obtained by the S2 inverse solution. 
Fig. 5. The blading of a multistage turbopump 










V = 0.7 

Fig. 6. Centrifugal impellers designed with r) = 1 and rj < 1 
having the same level of total pressure gain. 
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